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Abstract 

AdS/CFT correspondence is now widely used for study of strongly coupled plasmas, 
such as produced in ultrarelativistic heavy ion collisions at RHIC. While properties of 
equilibrated plasma and small deviations from equilibrium are by now reasonably well 
understood, its initial formation and thermal equilibration is much more challenging 
issue which remains to be studied. In the dual gravity language, these problems are 
related to formation of bulk black holes, and trapped surfaces we study in this work is a 
way to estimate the properties (temperature and entropy) of such black hole. Extending 
the work by Gubser et al, we find numerically trapped surfaces for non-central collision 
of shock waves with different energies. We observe a critical impact parameter, beyond 
which the trapped surface does not exist: and we argue that there are experimental 
indications for similar critical impact parameter in real collisions. We also present a 
simple solvable example of shock wave collision: wall-on-wall collision. The applicability 
of this approach to heavy ion collision is critically discussed. 
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1 Introduction 



The Quark Gluon Plasma(QGP) produced in Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Lab is believed to be strongly coupled [1] as evidenced by its rapid 
equilibration, strong collective flows well reproduced by hydrodynamics, and strong jet 
quenching. Applications of AdS/CFT correspondence [2] to strongly coupled QGP has 
generated many fundamental results [7, 8, 9, 10] , for some further results see [11] for a 
recent review. However this progress so far has been mostly related either to equilibrium 
properties of the plasma, or to kinetics/hydrodynamics close to equilibrium. 

The challenging issues related with violent initial stage of the collisions, in which the 
QGP is formed and equilibrated, producing most of the entropy, are not yet understood. 
One of them worth mentioning at the beginning is the strikingly different views on equilibra- 
tion held in statistical mechanics on one hand and in AdS/CFT-based dual gravity, on the 
other. Statistical/kinetic approaches treat equilibration and entropy production as some 
gradual deformation of particle distributions, from some initial non-thermal state toward 
the final cqTiilibrated one. In the dual gravity setting the source of temperature and entropy 
are both attributed to the gravitational horizons. Those may or may not be produced in 
a collision: for example decreasing the collision energy or increasing the impact parameter 
one may reach a point at which no horizons are formed. This implies certain singularities, or 
a view that switching in equilibration is similar to a phase transition rather than a gradual 
deformation. 

Formation of a black hole in a collision, which is then is falling toward the AdS 
center, was first considered in [18], with a spherical black hole. Janik and Peschanski 
[19] have proposed asymptotic (late-time) solution, corresponding to rapidity-independent 
(Bjorken) flow, see [20, 21] for most recent advances along this direction. 

Grumiller and Romatschke [3] tried to describe the initial stage of high energy col- 
lisions, starting with the gravitational shock waves of certain type. In section 5 we will 
explore the formation of horizon in similar setup, but taking a different point of view: the 
image on the boundary must be due to the source in bulk. This will lead to a different and 
more consistent initial conditions, as well as subsequent evolution of matter. 

A pcrturbative treatment of the initial conditions is discussed by Albacctc, Kovchcgov 
and Taliotis [5] . Other models of equilibration based on solutions to dynamical Einstein eqns 
include our previous work [22] in which a gravitationally collapsing shell of matter in AdSs 
space is considered. It sheds light on how formation of isotropic and homogeneous plasma 
may proceed through a very specific "quasiequilibrium" stage. We calculated the spectral 
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densities and found they deviate from their thermal counterpart by general oscillations. 
Another interesting solution describing isotropization of plasma was proposed by Chesler 
and Yaffe [23] recently. 

The issue wc will address in this work is formation of trapped surfaces and entropy 
production in the collision of two shock waves in AdS background. The work in this direction 
in AdS/CFT context had started with the paper by Gubser,Pufu and Yarom [4], who 
considered central collisions of bulk pointlike black holes. They had located the (marginally) 
trapped surface at the collision moment . Its area was then used as an estimate (a lower 
bound) of the entropy produced in heavy ion collision. In the limit of very large collision 
energy E they found that the entropy grows as E'^/^. In section 6 we will critically discuss 
how realistic are these results. 

In this work we extend their work in two directions. One is the obvious extension to 
collision of shock waves with nonzero impact parameter. We find interesting critical phe- 
nomenon, analogous to shock wave collision in Minkowski background [12, 13, 14]: beyond 
certain impact parameter, the trapped surface disappears and black hole formation does 
not happen. The other direction deals with a much simpler case of wall-on-wall collision, 
which was in a way overlooked before. 

2 Shock Waves Collision and Trapped Surface 

It is useful to review the main steps of [4] first. The AdS background can be written as: 



where u = t—x^ and v = t+x^. x^ is longitudinal coordinate and x^,x'^ are transverse 
coordinates. 

The shock wave moving in +x^ direction is given by: 



ds'^ = L 



2 -dudv + {dx^)'^ + {dx'^f + dz^ 
^2 



(1) 



ds 




+ L 



(^{x^,x^,z) 



z 




(2) 



with ^{x^,x'^ 



z) satisfies the following equation: 




(3) 
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The 5-Dimensional source Juu c&n be arbitrary function in principle. CH is the Lapla- 
cian in the hyperboUc space H^: 

,,^„, = L^i^^l±Spl±^ (4) 

The shock wave moving in —x^ direction can be obtained by the substitution u ^ v 
to (2) and (3). 

The marginahy trapped surface is found from the condition of vanishing of expansion 
^[15]. The trapped surface is made up of two pieces: 5 = 5i U 52- Si{S2) is associated 
with shock wave moving in +x^(— x^) direction before cohision. An additional condition is 
that the outer null normal to Si and 52 must be continuous at the intersection C = iSi n (S2 
point u = V = Q to avoid delta function in the expansion. 

To find out »Si associated with the first shock wave, 



as = L 2 o{u)au (5) 

the following coordinate transformation is made to eliminate the discontinuity in 
geodesies: 



v^v + -^Q{u) (6) 
where 6(n) is the Heaviside step function. Si is parametrized by: 



u = 0, V = —il>i{x^ jx"^ , z) (7) 

The expansion is defined by 9 = /i^'^V^/jy, with the outer null normal to Si. h^'^ 
is the induced metric. It can be constructed from three spacelike unit vectors w'^, W2, wt^, 
which are normal to Si: 

h'''' = w':^w\ + w^wl + w'^wl (8) 
The vanishing of expansion gives the equation: 



□ -^)(^i-$i) = (9) 
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with '^i{x^,x'^,z) = ^■ipi{x^,x'^,z). 

The vanishing of expansion on S2 associated with the second shock wave can be 
worked out similarly: 



(a - ^ (*2 - $2) = (10) 

At the intersection C = Si n (S2, Si and S2 coincide, therefore i{x^ ^ x"^ , z) = 
'^2{x^,x^,z) = 0. The continuity of outer null normal can be guaranteed by V^'i • V*2 = 4. 

In summary, the aim of finding marginally trapped surface becomes the following 
unusual boundary value problem: 



□ - ) (*i - ^1) = 



^-Jl) (^2-$2)=0 
*l|c = *2|C = 



(11) 



The boundary C should be chosen to satisfy the constraint: 



V*i-V*2|c=4 (12) 

Note (11) and (12) are written in the form of scalar equation, invariant under coordi- 
nate transformation. For central collision, the source Juu are identical for two shock waves. 
In [4], they are chosen to be 



Juu = E6{u)5{z - L)5{x^)5ix^) (13) 

The solution of $ corresponds to this source give rises to the following stress tensor 
on the boundary field theory: 



_ L2 ^(x^^x^,z)6{u) _ 2L^E 

" i3 - 7,(^2 + (^1)2 + (2,2)2)3'^ W (14) 

The special source (13) preserves an 0(3) symmetry in H^, which is manifest in the 
following coordinate system: 
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with the point source sitting at r = 0. We wih elaborate the symmetry later in the 
context of non-central collision. 

The 0(3) symmetry helps to solve (11) analytically. The area of the trapped surface 
can be calculated and give a lower bound to the entropy produced in the collision of shock 
wave, assuming the area theorem holds in AdS background. 

For non-central collision, the situation is complicated by the loss of 0(3) symmetry. 
In Minkowski background, the problem of non-central collision of point shock waves in 
D = A was solved beautifully in [12] by conformal transformation. In D > 4, it was solved 
numerically in [13]. In all critical impact parameter was found, beyond which the 

trapped surface seized to exist. 

In the next section, we will cast (11) into an integral equation, which allows us to 
solve (11) numerically. 

3 Calculation of the Trapped Surface 

Note (11) resembles the electrostatic problem in flat space, with being the electric poten- 
tial. We are familiar with the fact that the electric potential can be expressed as an integral 
of surface charge density. We want to see if this can be achieved in AdS space. 

Let us start with the electrostatic problem in flat space. Consider the following 
electrostatic problem, which is similar to (11): 

V^^iix) = V^^i{x) (16) 
^i{x)\c = (17) 
V^'i • V^2 = 4 (18) 

where i = 1,2, is the Laplacian in fiat space. *i is the electric potential corre- 
sponding to the source V^$i, placed inside an empty chamber with conducting boundary 
C. The boundary should be chosen properly such that the constraint (18) is also satisfied. 

We want to express the electric potential by an integral of the surface charge density. 
This can be done with the help of the free boundary Green's function defined as the solution 
to: 
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(19) 



with the solution given by: 

Take (16) multiphed by G{x,x') minus (19) multiplied by *i(x), and then integrate 
over the space inside C, we obtain: 



Denote Bi{x) = -2!™ (the magnitude of electric field on the boundary) and note 
^'i(x) vanishes on the boundary. With x' taken on the boundary C, (22) evaluates to: 



The constraint (18) is simply Bi{x)B2{x) = 4. Wc have cast a problem in the volume 
into a problem on its boundary C. (23) is a Fredholm integral equation of the first kind. 
We can use the following method to solve (16): Starting with some trial shape of C, we can 
solve (22) to obtain Bi{x) and check if (18) is satisfied. We can use iteration to tune the 
trial shape until (18) is satisfied. 

Now we hope to apply similar method to the problem of trapped surface, the difference 
being the space is instead of flat. 

As in case of electrostatic problem, we will keep using Green's function in AdS, deflned 
as the solution to the following: 



where g is the metric of H^. 

The Green's function was solved in [16, 17]. We quote the result here with L depen- 
dence restored. 





(23) 




(24) 
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1 e^" 

G{x,x')-- 



AttL sinh u 

cosh n = 1 + ' \ ^ (25) 



where u is the invariant distance in K^i^AdS'^). 
It also proves useful to note another relation: 



M 



= f -^d^{^g>'''d,f)^^e,pxdx''AdxfAdx^ 

= / di^g^^'^d.fe^px^dx'^Adx^) (26) 

where dx^ = g^^ yfge^px^dx'' Adx"^ . M is taken to be the manifold in bounded by 
C, the metric g refers to H^. f is arbitrary function of x. 

With (26) and (25) at hand, we are ready to proceed: 



(□ - ^) *,(x) = (□ - ^) $i(x) 
^ {n-^)Gix,x') = ^S('\x-x') 
with z = 1, 2. All the derivatives are with respect to x. The first line of (27) multiplied 
by G{x,x') minus the second line of (27) multiplied by ^'i(x), then integrate over M, we 
obtain: 

(g{x,x') (n - -^^ ^i{x) - (n - -^j G{x,x')^ ^d^x = 

G{x, x') (a - <^i{x)^d^x - ^iix') (28) 



J (g{x, x')d {dr,^i{x)dx''^ - ^i{x)d {d^G{x, x')dx''^^ 

G{x, x') (u - ^ ^i{x)^d^x - Mx') (29) 



J (d (g{x, x')d^^iix)dx''^ - d (^^i{x)d^G{x, x')dx''^^ 

G{x, x') {u - ^ <^i{x)^d^x - ^i{x') (30) 
/ (g{x, x')d^^i{x)d^'' - ^i{x)d„G{x, x')d^'') = 

JdM ^ ' 

G{x, x') {u - ^ ^^{x)^d'x - ^i{x') (31) 
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where in the last hne we have used Stokes theorem on manifold M. 
Putting x' on C, we can simplify the above with ^i\c = 0: 

/ Gix,x')d,^iix)d^" = [ G{x,x')(n-^]<^iix)^(fx (32) 

JdM JM \ ^ J 

Furthermore, we have d^il^dx" = on C since = 0. On the other hand, nvdx"\c = 
0, where rii^ is the unit vector normal to the boundary C. Therefore, we may write: 

d^^i = -BiU^ (33) 
With the help of (33), (32) and (12) can be further simplified to: 

- / G{x,x')Bi{x)dS= f G{x,x')(n-^)^i{x)^d^x (34) 

JdM JM \ ^ J 

Bi{x)B2{x)=4 (35) 
where dS = n^dx^ is the area element. 

Before proceeding to non-central collision, we would like to reproduce the 5-D result 
of [4] first. Working in spherical coordinates (15), the shape of C is parametrized by r = 
po = const. The simplest point shock wave corresponding to Juu = E6{u)d{z — L)6{x^)6{x^) 
is given by: 



4G,E 1 + 2(r/L)2 - 2r/LVl + (r/L)^ 
= $2 = — ^ (36) 

The Green's function (25) is invariant under coordinate transformation. In spherical 
coordinate, it is given by: 



1 e^" 

Gix,x')-- 



AttL sinh u 

cosh u = ^J7^Jl?T\^J7^liy+\ - rr'/L^ (cos 9 cos 6' + sin sin 6' cos((/> - 0'))(37) 

In the presence of 0(3) symmetry, it is sufficient to show (34) holds for 0' = 0, when 
the integral in cj) is trivial. On the other hand, (35) implies B\ = B2 = 2. As a result, we 
only need to verify: 



27r r de{-2) ^ -/"^ pI sin 9 = ^^^^ " (-4^5^)4^ 
Jo smhn pQ 

cosh u = Po + 1 ~ Po cos 9 (38) 
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It is not difficult to complete the integral in 0, we finally arrive at 2G5E = -^/TTp^TZ^Po' 
which is equivalent to (115) in [4]. 

4 Colliding Point Shock Waves at nonzero Impact Parameter 

4.1 Shock Waves in Spherical Coordinate 

Consider two shock waves with impact parameter b, given by: 



(39) 



The corresponding stress energy tensor associated with two shock waves are given by: 



"""7r(z2 + (xl-|)2 + (x2)2)3'^W 



2L^E 
niz^ + (xi + 1)2 + 

Therefore zq characterizes the size of the nucleus. Wc will use spherical coordinates 
in solving (34). In case of central collision, when 6 = 0. The shock wave center can be 
placed at the origin of spherical coordinates r = 0. This is achieved by first going to global 
coordinates ¥"^{1 = 0, 1, 2, 3): 



2 



2 V z 



Y^ = L- 



2 V " ' 
1 



= L— (40) 



z 



The global coordinates link to spherical coordinates in the following way: 



yo = + L2 
= r cos 9 

Y"^ = r sin 9 cos (f> 

= r sin 6* sin ^ (41) 
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When 6 = 0, the center of the shock waves can be put at the origin if we set A; = 
The possibiHty of moving any point to the origin reflects the maximally symmetric property 
of AdS space. 

When 6 7^ 0, wc want to place the two shock waves at opposite positions with respect 
to the origin, so that the boundary of trapped surface C will have axial symmetry. Setting 
1 + ^ = -p^, we have Y"^ = = and Y^ = According to (41), we have the shock 

waves at r = ^ = and 9 = tt. The differential equation in (34) becomes: 

(n -^)^i = -16nG,E^^^^£^S{r - ro)<5(0 - 9,)Sicf>) (42) 

where vq = i^, 9i = 0, 02 = t^- We observe that in spherical coordinate, the trapped 
surface only depends on G^E^ and tq. Since AdS radius is a free parameter, which will 
not appear alone in the final result in dual field theory, we may set zq = L without loss of 
generality. As a result we have h = 2ro. 

4.2 More General Shock Waves 

Before proceeding to numerical study of trapped surface, we choose to take a moment to 
investigate the symmetries of the problem, which will help us to study more general shock 
waves. To see this, we prefer to work in the differential form of the problem: (11) and (12). 

As we noticed before, (11) and (12) are scalar equations. *i is a scalar. It is invariant 
under coordinate transformations: x ^ x, ^i{x) ^i{x) the boundary remain the same 
C ^ C, but takes a different functional form in new coordinate. As a result the third line 
of (11) and (12) are automatically satisfied. Suppose the transformation also preserves the 
form of the operator: □ — then ^i{x) becomes another solution to (11) and (12). 

We will focus on transformations that leaves the center of the shock waves on the axis of 

6 = 0, TT. 

To identify such a coordinate transformation, we first make a change of variable: 



rsm9 = t 

r cos 9 = a/L^ + sinh rj 



The metric of becomes: 



d^' = + + t')dri' + t'dcj^' (43) 

1 + r^/i^^ 
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The metric is 77 independent, therefore the transformation: t = t,(f) = (f),fi = 77 + Arj 
will not change the operator □ — 3/L^. i = t also guarantees the center of the shock waves 
remain on the axis of 9 = 0, tt. We have obtained the desired coordinate transformation, 
which is just a translation in t]. It is easy to work out the corresponding transformation in 
spherical coordinate: 



f sin ^ = r sin ^ = f 

f cos 9 = \/ly^+t^ sinh(?7 — Arj) 

r cos 9 = \jL?^r^ sinh 77 (44) 

One can verify explicitly (44) preserves the form of (15). (44) moves the center of 
the shock waves from _ yS _ yl _ y2 _ y3 _ yl _ _|_^^ pQgj^ _ 

a/Z^ + Tq sinh Ary. This means collision of shock waves centered at Y'^ = = 0, = 
ibro cosh At] — Y^L^ + Tq sinh At] will generate the same entropy as those centered at Y^ = 
Y"^ = 0, Y^ = zbrQ. This allows us to study the collision of more general shock waves. Let 
us consider the following shock waves: 



(n - = -167rG5Eu6iu)6{z - Zu)S{x^ - Xu)S{x^) 

(n - ^2 = -16TrG5EJ{v)S{z - z^)S{x^ - x^)S{x'^) (45) 

In this paper, wc restrict our interest in shock waves with identical invariant energy, 
defined by = E^^ = E. This keeps the mirror symmetry of the problem intact. We 

will see the center of the shock waves can be placed at = Y^ = 0, Y^ = ±ro cosh Arj — 
\/L2 + sinh A?7. This is equivalent to the statement that a solution can always be found 
to the following equations: 



L— = tq cosh A?7 — -i/ + rn sinh Ar? 
Zu V " 

L— = — ro cosh A77 — \/ + r?. sinh A77 

Zv V " 



L 



Ax 



(46) 
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(46) can be solved easily by switching to the variable rjo = sinh' 



^ ^. A solution to 



(46) always exists for any given Zu, Zy and Aa;. We include the corresponding ro here, as it 
is the only relevant quantity, apart from E, for entropy calculation: 



In summary, we have shown that for shock waves (45) satisfying = Ey-^ = E, 

the entropy is only a function of G^E and (47). Note ro is only a function of invariant 
distance between the center of shock waves. 

4.3 Numerical Solution of Trapped Surface 

With all the simplification, we are ready to find the trapped surface for different impact 
parameter. Our procedure is as follows: Axial symmetry allows us to parametrize C by 
r = p{0). Integral in <p on the LHS of (34) can be expressed in terms of elliptic integrals. (34) 
becomes essentially 1-D integral equation. We discretize the integral by 199 points, equally 
spaced in the full range of 6. The integral on the LHS of (34) is discretized accordingly, 
and the integral on the RHS can be expressed in terms of elementary function due to the 
simple form of the shock wave. We use the same sample points for 9', bringing (34) into a 
matrix form: 



where the indices i,j = 1, • • • , 199. K{6j,9'^) contains the Green's function and the 
induced metric. S{6'j) is from the RHS integral of shock wave. 

A special treatment is needed for diagonal matrix element of K{9j, 9'^) where 9j = 9[. 
The explicit integrand expressed in terms of elliptic integrals shows the it is logarithmically 
divergent in \9 — 9'\, yet the integral is convergent. The integral in this interval, represented 
by the diagonal matrix element, is estimated by sampling the integrand by certain number 
of points in the interval. The sample integrand are used to extract the coefficients of terms 
\n\9 — 9'\, 1, {9 — 9')\n\9 — 9'\ and 9 — 9' hy method of least squares. Those coefficients are 
finally used for calculation of diagonal matrix elements. 

The mirror symmetry of the two shock waves implies D2{0) = Bi{-k — 9). There- 
fore it is sufficient to calculate one of them. Given a trial shape of trapped surface 




(47) 




(48) 
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Tabic 1: critical impact parameter at different energies 
^ 0.1 0.5 1 4 9 12 15 50 100 

^ 0.40 0.86 1.14 1.90 2.50 2.74 2.94 4.28 5.30 

r = p{9), which is also necessarily symmetric under 6 ^ n — 6, we can solve for B(9) 
from (48). We then evaluate A(6') = Bi{9)B2{9) — 4 and tune the shape function ac- 
cordingly. We repeat the process until (35) is satisfied to certain accuracy. In order to 
assure fast convergence, we find it very helpful to calculate the gradient of p{9). The 



gradient is the matrix form of the functional derivative: ^^gj^- We parametrized C by: 
p{9) = On cos 2(n — 1)9, where M is a truncation number. The same decomposition 

applies to A(^): A(^) = Yln=i ^" ^('^ ~ ^)^- '^^^ gradient in this representation is given 
by a M X M matrix: which again contains elliptic integrals. For a given collision 

energy, we can find the boundary C until certain critical impact parameter is reached. The 
critical impact parameter is located where ^^p- diverges[13]. Empirically, the gradient 
converges as A{9) reduces in the iteration, if the impact parameter is within the critical 
value. The gradient diverges as A(^) reduces in the iteration, if the impact parameter lie 
beyond the critical value. 

Fig.l shows the shapes of trapped surface at = 1 and = 100 for different 
impact parameters. The shapes are represented in spherical coordinate. We observe the 
critical trapped surface does not scale with collision energy in spherical coordinate. As 
collision energy grows, the trapped surface gets elongated in the axis of mismatch and 
larger M is needed to reach prescribed accuracy. 

We also obtained several critical impact parameters corresponding to different ener- 
gies. The results are listed in Tabid 

Fig. 2 shows the log- log plot of critical impact parameter versus collision energy. It 
suggests a simple power law within the energy range used in the numerical study. The data 
are fitted with ^ = a (^)^ to give: 



a = 1.07, p = 0.37 (49) 

b ~ E^L^~'^/^, the numerical value from fitting shows the critical impact parameter 
grows with collision energy and nucleus size. 

The area of the trapped surface(twice the area of C) sets a lower bound of the entropy 
produced, given as follows: 
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theta- 




Figure 1: (left)The shapes of C (the trapped surface dX u = v = 0) at ^ = 1. The impact 
parameters used in the plot are 0.4L, 0.6L, 0.8L, l.OL, I.IL, 1.14L from the outer to the 
inner. The innermost shape being the critical trapped surface. (right)The shapes of C (the 
trapped surface at n = u = 0) at = 100. The impact parameters used in the plot are 
l.OL, 2.0-L, 3.0L, 4.0L, 5.0-L, 5.3-L from the outer to the inner. The innermost shape being 
the critical trapped surface. As collision energy grows, the trapped surface gets elongated 
in the axis of mismatch. 



2A 1 

'trapped = = 7^ / y/Q^^X (50) 



2G5 

where A is the area of the boundary C. The prefactor is 



- (51) 

G5 TT 

We plot the lower bound of entropy in the dual field theory for energy = 100 in 

Fig.3 



5 Wall-on-wall collisions 

In this section we address a simpler form of the shock waves, called wall-on- wall in [18], in 
which there is no dependence on two transverse coordinates. Grumiller and Romatschke [3] 
have also discussed it, using gravitational shock waves. The problem with their approach 
is (power) growing amplitude of the shock as a function of holographic coordinate z. If 
so, collision dynamics resembles the atmospheric turbulence in the sense that the largest 
perturbation is at the largest z - namely in the infrared modes - cascading down toward 
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Figure 2: The log- log plot of critical impact parameter versus collision energy. 



Area/L'^3 




: : : : , : : : : , : : : : , : : : : , : [j/L 

2 3 4 5 

Figure 3: The scaled entropy 2^56'/^^ (the area of C) as a function of the impact parameter 
scaled b/L. The energy used is = 100, where ^ = 

higher momenta (UV). First of all this is physically different from the heavy ion collisions, in 
which case the initial wave function have partons well localized near the so called "saturation 
scale", from which the equilibration domain (trapped surface) propagates both into small 
z (UV) and large z (IR) directions. Second, we think it is also inconsistent: a function 
growing with z cannot be considered as a small perturbation to the background metric, 
which is decreasing at large z as l/z^. One may think that gravity near the AdS center 
z = 00 should never be touched, as this is the original positions of the branes which 
gave the basis to AdS/CFT correspondence in the first place. 

Our choice of the initial conditions, describing colliding walls with fixed parton density 
and thus fixed saturation scale is given by the following source 

(a - = -Wn^6{z - zo) (52) 

The corresponding solution to Einstein eqn subject to the boundary condition "i^iz) ^ as 
z — > is easily obtained: 
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^0 



(53) 



Note that it decreases in both direction from the original scale zq: therefore (as we will see 
shortly) the trapped surface has finite extensions in both directions from it. 

The corresponding stress energy tensor on the boundary (as seen by an observer living 
in dual gauge theory) is 

EL'' 



T — 

11.11. — 



-5{u) 



(54) 



The stress energy tensor is the same as that used in [3, 5], and our solution converges well 
as 2: ^ 0. We choose to collide states with different energy therefore we fix zq^ but use 
different E. Applying the general discussion of shock wave in Sec. 2 and noting the trapped 
surface only depends on z, we obtain: 



zH'l - z-^'i - 3^i = -lQTTG^Ei5{z - zq) 

^i{Za) = ^i{Zb) = 



(55) 
(56) 
(57) 



with i = 1,2. ^'i and ^'2 arc shape functions corresponding to two shock waves. The 
trapped region at u = v = Q is limited by the interval Za < z < zjj. The constraint (35) 
takes a simple form: 



*l(^6)^'2(^(,)^=4 



(55) is easily solved to give: 



(58) 



D 



C = -AirGsEi 



C{{^] -^j z<zo 

, 3 



(i) - 



Z > Zq 



(s) - (t) 
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\3 

D = -AirG^Ei >^ '-^^^ (59) 



Plugging (59) in (58), we obtain: 



87rG5/grg2 .... 

2a + ^6 = ^ (60) 

Note £^1, £^2 appear only in the combination \/E^E2, This is consistent with the 
picture that only the center of mass contributes to the entropy. Recall the the center of 



mass of two massless particles with energy Ei, E2 is 2\/EiE^. The resulting cubic eqn (60) 
can be solved by Cardano formula, but the explicit solution is not illustrative and is not 
showed here. The entropy is given by: 



^ 2A_ ^ J ^dzd^x± 
S ,1 1 . 

The leading behavior of entropy per transverse area s in energy is extracted: 



s ^{7rG5VEiE2zi)3 (63) 

^0 

The power 2/3 is the same as point shock wave obtained in [4]. There is also an 
obvious lower bound of the energy for the formation of trapped surface: 



AnG^E > zq (64) 

The equality is reached at Za = zi,, when the C has vanishing volume. For general 
energy, s is evaluated as a function of effective colliding energy E = -v/E^Ej- We again set 
Zq = L. Fig.4 shows the entropy as a function of effective colliding energy. 



6 Matching heavy ion collisions to those of gravitational shock 
waves 

As the reader who came to this point of the paper knows, it was up to this point a methodical 
work devoted to solving well-posed mathematical problems set in the AdS/CFT dual-gravity 
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Figure 4: The scaled entropy per transverse area ^ (the area of C per transverse area) as 
a function of scaled effective colliding energy G^E jl? ^ where ^ = 

framework. However now, near the end of this work, we would like to address a wider issues 
of applicability limits of such approach, as well as the best strategy to use it for practical 
problems. 

Gubser et al [4] have applied the gravitational collision scenario literally, selecting 
initial conditions at time long before nuclei collide. More specifically, they have (i) tuned 
the scale L or zq of the bulk colliding object to the size of the nucleus R and (ii) have used 
the realistic CM gamma factor of the colliding nuclei E/m = 7 ~ 100. The result of such 
choice is a completely unrealistic fireball produced, in spite of a reasonable entropy. Indeed, 
the size of the trapped surface [4] is huge, about 300/m, which is very large compared to 
colliding nuclei. In real heavy ion collisions the produced fireball has the same size as the 
nuclei, with the radius about 6fm. The initial temperature — as estimated by Zmin ~ 
where Zmin is the minimal distance of the trapped surface to the AdS boundary - is however 
way too high. So, what went wrong with this straightforward approach? 

The answer to this question is in fact well known: initial formation of the partonic 
wave function, describing nuclei at the collision moment, can not be adequately described by 
the dual gravity. We know from experiment that growing partonic density makes hadrons 
and nuclei blacker and of larger size, as the collision energy grows. This is usually described 
by a Pomeron fit in which cross section ~ Although qualitatively similar to what 

happens in gravitational collisions, this growth is very compared to that predicted in dual 
gravity. Indeed, the observed Pomeron intercept a{t = 0) — 1 ~ 0.1 while in the AdS/CFT 
world the Pomeron intercept a{t = 0) — 1 = 1[24] . Thus the effective size of objects in 
gravitational collision grows with energy with an exponent ten times that in the real QCD. In 
view of this, one should clearly give up the idea to tune the scale L or zq of the bulk colliding 
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object to the size of the nucleus, and tune it perhaps to the parton density ("saturation 
scale" Qs in the "color glass" models) of the corresponding nuclear wave function. 

More generally, we are dealing with a complicated problem in QCD, in which the 
effective coupling runs, from higher scale to lower as the collisions progress from initial 
violent partonie stage toward equilibration, expansion and cooling. So in principle, it would 
be logic to switch - as smoothly as possible - from the weak-coupling based methods (such 
as classical Yang-Mills) to strong coupling ones (such as AdS/CFT) at certain proper time 
Tgwitch appropriately chosen by the evolution of the coupling*^. 

Therefore one should not try to tune the parameters of the gravitational collision 
model neither to initial nuclei, at r = — oo, nor to "decoherent" partons at the collision 
moment, at r = 0, but at the later time TgwUch-'^'^ Although we at the moment do not 
understand the evolution of appropriate coupling quantitatively enough, one may always 
treat it as a parameter. The practical utility of the AdS/CFT approach at later time r > 
Tgwitch still remains significant: namely one can use much more fundamental dual gravity 
description instead of its near-equilibrium approximation, the hydrodynamics, currently 
used. 

7 Are there critical impact parameters in heavy ion colli- 
sions? 

Summarizing our findings in one sentence, those is the existence of the discontinuity in 
grazing gravitational collisions in the AdS space. As one smoothly increases the impact 
parameter b, the trapped surface and black hole formation disappear suddenly, at certain 
critical impact parameter bc{E) depending the collision energy E. The reason for this seems 
quite general: increasing b one increases the angular momentum of the system while at the 
same time decrease the mass which can be stopped, and at some moment — as one knows 
from Kerr solution for rotating black holes ~ black hole formation becomes impossible. 

Suppose the AdS gravitational shock waves can describe the strongly coupled plasma 
in heavy ion collisions: then one would expect similar behavior in heavy ion experiments. 
We have looked at the data and found that indeed there are experimental indications that 
relatively rapid switch of the underlying dynamics at some bc{E) seem to exist. 

*'^The so called AdS/QCD approach (see e.g.[31, 30]) tries to incorporate the running coupling into the 
gravitational framework. A particularly simple example of that is a jump of the coupling at certain "domain 
wall" scale proposed in [33]. 

*^It is proposed in [5] that one may choose to collide some special unphysical shock waves 
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Figure 5: (left) PHOBOS data [25] on integrated number of charged particles, scaled by 
Npart/'^^ in P+P) d+Au, Cu+Cu and Au+Au collisions as a function of centrality. The 
uncertainty of Npart has been included in the error bars, (right) The height of the ridge as 
a function of the number of mean binary collisions per nucleon. The data are from STAR 
collaboration [27, 28] at two collision energies shown in the figure. 

The most straightforward observable is entropy, related to the particle multiplicity 
versus the impact parameter. In Fig.5(left) from [25] we show some data plotted as a func- 
tion of the number of participants Npart- The right end of the figure corresponds to all 
nucleons participating, or central collisions: toward the left end are peripheral collisions. 
There are indeed two values of multiplicity per participant observed, one for small systems, 
pp and dAu collisions (stars and crosses) and one for "large" systems, CuCu and AuAu 
(circles and squares). There must be a transition between them somewhere, but, unfortu- 
nately, the experimental multiplicity measurements for "grazing" collisions are not available 
yet So, unfortunately, we do not yet know how exactly transition from one regime to 
another happens and what is bc{E), if it can be defined. 

However some other observables associated with collective fiows of excited matter do 
show rapid changes at certain hc{E) seems to be there. Some evidence for that were seen 
in the elliptic fiow measurements, as deviations from the hydro dynamical predictions for 
very peripheral collisions. Even more clearly those are seen in the centrality dependence 
of the so called "ridge" phenomenon (see its relation to fiow in [26, 28]) which we show in 
Fig.5(right). 

Admittedly, these rapid change of the dynamics have not been systematically studied 

*^Small multiplicity collisions are detected for all systems, but their accurate separation from beam-residual 
gas collisions has not yet been systematically resolved. 
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yet, neither experimentally nor theoretically. The naive explanation often given to it at- 
tribute the change to the fact that it happens when overlap system gets "too small" in terms 
of participating nucleons Np, causing large enough fluctuations 0(1/^^^). However, if this 

would be the reason, one would expect this jump to be dependent on Np and independent 
on the collision energy. Furthermore, the gravitational collisions do not have any discrete 
elements at all, while predicting hc{E) growing with E, as observed in Fig. 5. We therefore 
suggest that angular momentum may also be important: this issue clearly deserves to be 
studied further. 

8 Conclusions 

In this work we have developed a method to solve for the shape of the trapped surface 
based on an analogy to electrostatic problem in flat space: its main idea is to proceed from 
differential to integral form of the equation. We used the method to obtain the shape of 
trapped surface at different impact parameters and collision energy. We observe a critical 
impact parameter within the range of energy we explored. The phenomenon is analogous 
to the the critical behavior found in flat space[12, 13, 14], the difference being the critical 
trapped surface depends both on the collision energy and the nucleus size. We found the 
dependence is approximately given by a power law. Furthermore, the shape of the critical 
trapped surface gets elongated in spherical coordinate as the collision energy grows. We 
also discussed in the preceding subsection that grazing heavy ion collisions also seem to 
suggest a rapid switch to another dynamics, without equilibration. The exact cause of this 
jump is to be clarified in further studies. 

We also studies wall-to-wall collision of shock waves as a simple version of the problem. 
The wall is sourced by a delta function at certain initial scale Zq. We believe it is more 
reasonable initial conditions than those used by Grumiller and Romatschke [3] , to be used 
in future following their method to study the initial stage. 

The applicability and limitation of this approach is discussed. We claim it is more 
realistic to adopt partonic picture in initial stage and only switch to effective gravity treat- 
ment at some time after collision, when the coupling becomes strong enough. However, we 
argue that the observed critical phenomenon is still relevant for heavy ion collisions, where 
there also seems to be rapid change of collision regime as a function of impact parameter. 

Finally we would like to mention very recent work by Alvarez-Gaume et al[29] who 
discussed another extension of the problem. They considered central collision of shock 
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waves sourced by certain nontrivial matter distribution in the transverse space. They in 
particularly discuss critical phenomenon occurring as the shock wave reaches some diluteness 
limit and the formation of the trapped surface is no longer possible. It would obviously be 
interesting to study how the two forms of critical phenomenon arc related. 

Note added in version 2. Horatiu Nastase had informed us about his early work [35] 
in which he also addressed non-central gravitational collisions, see also summary in [34]. 
His estimate of the critical impact parameter bc{E) E^/^ at large E is quite different 
from our power, pcrhapsc because we are not working at asymptotically large E. We were 
also provided by (so far not posted) work by the Princeton group [36], who compared our 
numerical results for noncentral collisions with their analytical formulae and observed good 
agreement. 
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